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A recently introduced model of dually weighted planar graphs is solved in terms 
of an .elliptic parametrization for some particular collection of planar graphs de- 
scribing the 2D R 2 quantum gravity. Along with the cosmological constant A one 
has a coupling f3 in the model corresponding to the ii 2 -coupling constant. It is 
shown that for any value of (3 the large scale behavior of the model corresponds to 
that of the standard pure 2D quantum gravity. On small distances it describes the 
dynamics of point-like curvature defects introduced into the fiat 2D space. The 
scaling function in the vicinity of almost flat metric is obtained. The major steps 
of the exact solution are given. 

1 Introduction 

Matrix models counting the number of planar graphs of various types have 
played an important role in various domains of physics and mathematics, es- 
pecially, in the quantitative approaches to two-dimensional quantum gravity 
and non-critical strings and string field theories. 

By means of a standard NxN hermitean 1-matrix with the partition func- 
tion: 



we generate (by expanding logZ in powers of tj.) all the abstract graphs 
weighted with the product of factors tk each corresponding to a vertex with 
the coordination number k and with the overall factor N 2 ~ 29 where g is the 
genus of the graph. The t^s allow us to control the frequencies of coordination 
numbers of the vertices, but not the coordination numbers of the faces. One 
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where the potential is defined by 
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of our main technical goals is to generalize the 1-matrix model ( [l]) so that we 
can also weight the faces with independent weights. 

The model (Q) has been proven to describe (with the appropriate choices 
of ifc's) the multicritical points of the pure 2D gravity corresponding to the 
rational conformal matter (of the type.l2,2n-l) in the standard classification 
of the rational 2D conformal theories) BB. 

An interesting question to ask is whether the 1-matrix model could be 
used for the description of the pure 2D gravity with higher derivative terms, 
described by the following formal functional integral: 

Z{\, 0) = J Vg ab (z) e~ / d2z V^~9 (p+aJ?. 3 +^fl 2 9 +...) ^ 

On the first sight, from the simple dimensional analysis, apart from the 
cosmological term (controlled by fi) and the (topological) Einstein term pro- 
portional to the curvature R g , it does not seem meaningful to put further terms 
into the action of 2D gravity. The simplest term one might want to consider is 
-j^Rg- The bare coupling constant 0o is however dimensionful and thus should 
be proportional to the cutoff squared. So it is small and in principle should be 
dropped, as well as any further higher derivative terms. 

On the other hand, if we make smaller and smaller the characteristic 
metrics in the functional integral ( ||) should approach the flat one, since the 
R 2 suppresses the fluctuations of the metric. So, an interesting question is 
whether one can get some interesting non-perturbative behavior, where the 
R 2 coupling could play some essential role. In other words, the question is 
whether there could be a non-perturbative phase of almost flat metrics, beyond 
some hypothetic "flattening" phase transition. It is clear that all the known 
approaches (like the Liouville theory approach) which start from the continuum 
formulation ( ^ are not valid for this purpose since the formulation might need 
serious modifications in this non-perturbative regime. Some nonperturbative, 
lattice formulation is thus needed. 

Unfortunately, the model ( ^), in spite of the infinite number of indepen- 
dent couplings, cannot be used to study a flattening transition. It is impossible 
to tune the couplings in such a way that we only generate flat, regular lattices 
corresponding to a flat metric. The model ( [j]) provides no control over the 
occurence of different types of faces, so even for only one non-zero coupling, 
say, tf. = <Jfc 4 5 we get a sum over arbitrary 4 graphs, which describes the 
highly developed fluctuations of pure 2d quantum gravity. 

Hence we have to introduce a model where we can control the numbers of 
faces with a given coordination number. The most general model of this type, 
containing the second infinite set of dual couplings which weight the faces 
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(dual vertices) of graphs, is given by the following partition function (here we 
fix the genus of the connected graphs g = 0): 



z(*v)=z n tf v; t* v « (4) 

G V*,v 9 eG 

where v* , v q are the vertices with q neighbours on the original and dual 
graph, respectively, and #i> q are the numbers of such vertices in the given 
graph G. This expansion is generated by the following matrix model: 

Z(t*,t) = J VM e~% tr M2 + Ntr V *( MA ), (5) 

with 

Vb {MA) =Y\ \ trB k (MA) k . (6) 

k=l 

A and B are fixed external matrices encoding the coupling constants through 

r q = \jftr* and t q = i I tr A*. (7) 

The model generalizes, for A ^ 1, the standard one matrix model ( [l]). Wc 
will call it the model of dually weighted graphs (DWG). 

The DWG model provides the possibility of flattening of the typical graphs. 
If we take, say, i& = t* k = S^a we will single out the regular square lattice 
in the partition function. The choice = 6^.3, t* k — <5fc,6 brings us to the 
regular hexagonal lattice, dual to the regular triangular lattice. These limits 
correspond to the (3 = limit in the continuum action of ( ||). Our purpose 
here is to investigate the critical behaviour of the DWG-model as we approach 
closer and closer this limit. 



2 A Solvable Model of R 2 2D Gravity 

In the next section we will show that the DWG-model is in principle solvable 
for any choice of coupling constants {t, i*}. But for our physical purposes we 
don't need so big space of independent couplings. We have to find the simplest 
DWG-model with only two independent couplings, one corresponding to the 
cosmological coupling A, another to the R 2 coupling (3. 

We have found that the simplest model of this type corresponds to the 
following choice of couplings: 

t 2 = VA t, U = A, t 6 - A* ^, ... t 2q = At {^) (q - 2) . (8) 
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With these, weights, it is easy to prove, using Euler's theorem, that the 
partition sum a becomes 

Z(t,\,/3) =i 4 XA !3 2( * V2 - 4 \ (9) 

G 

where A is the number of plaquets of the graph G and the number of pos- 
itive curvature defects. Note that the latter are balanced by a gas of negative 
curvature defects, whose individual probabilities are given in ( ||). 

After tuning the bare cosmological constant A (controlling the number of 
plaquets) to some critical value A c (/3), we expect this model to describe pure 
gravity in a large interval of (3. On the other hand, for A fixed and (3 = wc 
entirely suppress curvature defects except for the four positive defects needed 
to close the regular lattice into a sphere. It is thus clear that j3 is the precise 
lattice analog of the bare curvature coupling /3 in ( |||). The phase (3 = 
of "almost flat" lattices - very different from pure gravity - was discussed in 
detail in 

3 Sketch of Solution 

3.1 The Itzykson-DiFrancesco formula for the DWG-model 

The basic fact which allows to solve the DWG-model is the representation of 
the matrix integral ( |5|) in terms of the character expansion with respect to 
the irreducible-.representations R of the group GL(N), given by Itzykson and 
DiFrancesco tHl: 

_ TT (h e - lV'/i°" 

*(*.*•) = « £ n\e_; o V" x {h} (A) X W (B). (10) 

The sum runs here over all representations R{h} characterized by the 
Young tableaux with hi — i + 1 boxes in the i-th row. The non- negative integers 
hi obey the inequalities: 

h i+ i > hi. (11) 

Only the representations with equal number of even(odd) weights h e (h°) enter 
the sum in ( [l(]) . 

The Weyl-Schur characters are defined in the standard way: 

X{h} (A) = det(P hk+1 ^(9)) (12) 

(fe,o 

as a determinant of Schur polynomials P n {6) defined through 
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oo 

J2 z " p n{0) with 6i = -tr[A% (13) 



i 

n=0 



3.2 The Saddle Point Equation for the Most Probable Representation in the 
Planar Limit 

With the Itzykson-DiFrancesco formula we achieve what is necessary to make 
the large N limit analytically treatable: the drastic reduction of the number 
of degrees of freedom. Instead of N 2 integrations over the hermitean matrix 
Mij we are left with only N summations over the highest weight components 
hi < h% < ... < /ijv- If we assume that the the characteristic h's are of 
the order N (which will be confirmed by the solution) we can see that the 
summation factor in ( [To| ) is of the order exp N, and as usually one can apply 
in this situation the saddle point approximation. 

For our particular model ( ^- ^ the saddle point equation (SPE) is obtained 
by equating to zero the logarithmic derivative of the weight of summation in 
( |l^) with respect to h e or h° (we assume identical distributions for h e and 
h°). The SPE reads: 

2F(h) + P [ dh' -f^-j- = - In h. (14) 
Jo h h 

where 

Fth ^ 2 w t kJ K(br- < 15 > 

and p(h) is the density of the highest weight components in the represen- 
tation R. It is related to the resolvent H(h) of the highest weight components: 



t^h-hk .In h-h! 



m = Y.ir L ur= *'f^- (16) 



Although the details of the derivation of the SPE contain some subtleties, 
we can roughly say that the term 2F(h) comes from the two characters in 
( |l^) normalized by the Van-der-Monde determinant A(h) — Yli>j{hi ~~ 
(one takes A = B for the model dually equivalent to our model), the second 
term comes from log derivative of A(h) and of the product J\ i j(hf — h°), and 
log ft, comes from the Sterling asymptotics of (h e — 1)!! (or h°ll). 

As has been discussed in detail in □, this equation actually does not hold 
on the entire interval [0, a], but only on an interval [6, a] with < b < 1 < a: 
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Assuming the equation to hold on [0, a] would violate the implicit constraint 
p(h) < 1 following from the restriction ( [ll]). The density is in fact exactly 
saturated at its maximum value p(h) = 1 on the interval [0, b]. 

A similar saddle-point method for the representation expansion was first 
successfully used in for the calculation of the partition function of the mul- 
ticolour Young-Mills theory on the 2D sphere (see also B : u). 



3.3 Calculation of Characters in the Large N Limit 

To solve the SPD ( Q) for H(h) we have to find a method to calculate ef- 
fectively the characters, i.e. the function F(h). In this section we describe a 
method to do it, which is based on the identities for Schur characters, based 
on the basic properties of Schur polynomials. Leaving^ aside the details of the 
derivation (which can be found in our basic paper EJ), we summarize these 
identities in one equation already in the large N limit. Introducing the func- 
tion: 

G(h) = e H ^ +F ^ (17) 

we have: 

Q , CO „ o 

q=l 9=1 H 

where the coefficients of the positive powers of G in ( |l8|) are directly related 
to the correlators of the matrix model dual to the model defined by the eqs. 
( §- §), i.e. the model with potential V Ai {MA) = \ (MA) 4 : 

2, a L ^ ,i 



N dt 



2<l 



(x m {a))=(^tr{MAfi). (19) 



We have also assumed for the moment that only a finite number Q of 
couplings are nonzero (i.e. ti q = for q > Q). Furthermore, we were able to 
show in our work El that ( [l8| ) implies the functional equation 

e*W = LiL ±l[G q (h), (20) 

where the G q (h) are the first Q branches of the multivalued function G(h) 
defined through ( |lj) which map the point h — oo to G — 0. The saddlepoint 



equation (14), together with ( g0|), defines a well-posed Riemann-Hilbert prob- 



lem. It was solved exactly and in explicit detail in □ for the case Q — 2, where 
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the Riemann-Hilbert problem is succinctly written in the form 



2ReF(h)+H{h) = -ln(-A) 

2F(h) + ReHQi) = -ln/i, 1 ' 

where ReH(h) denotes the real part of H (h) on the cut [b, a] and ReF(h) 
denotes the real part of F(h) on a cut [—00, c] with c < b. This case corresponds 
to an ensemble of squares being able to meet in groups of four (i.e. flat points) or 
two (i.e. positive curvature points). We termed the resulting surfaces "almost 
flat" . It turned out that all the introduced functions could be found explicitly 
in terms of elliptic functions. 

The method can be easily generalized to any Q. However, the solution 
cannot be in general expressed in terms of some known special functions. In 
the next section we will show that the model of our interest can be solved in 
terms of elliptic parametrization, more complicated than its particular case, 
Q=2, of almost flat planar graphs. 



Let us also note that the equation ( 20) can be explicitly solved for F(h) in 



terms of a contour integral for some wide class of big Young tableaux. Hence 
we can find the character rather explicitly knowing H (h) and a set of Schur 
constants (t\, t<i, to) (see a). 

3.4 Solution of the Lattice Model of 2D R 2 Gravity 

Let us notice that for the model defined by ( |q, |S ) which we chose as a lattice 
realization of the 2D R 2 gravity, the equation ( |1S ) contains an infinite number 
of both negative and positive powers of the G-expansion. 

Labeling the first two weights as in u: t 2 , £4, we now include all the even 
t2 q with q > 3, assigning them the following weights: t2 q = t/±t q ~ 2 . Equation 
( |lj) can then be written compactly as: 

h-l = ^ + j£ r + positive powers of G. (22) 

G G (Cx — e) 

Dropping the positive powers of G (for G small enough) and inverting this 
equation, we see that there are two sheets connected together by a square root 
cut running between two finite cut points, d and c. On the physical sheet a 
further cut (running from b to a), corresponding to e H ^ h \ connects to further 
sheets. 

The equations ( |2l|) now read: 

2ReF(h)+H(h) = -ln(a^) ^ 
2F{h) + ReH(h) = -Ink, V ' 
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the first coming from the large N limit of the character ( ^) and the second, 
in view of ( |l6|), being the saddlepoint equation ( |l4|). These two equations 
tell us about the behaviour of the function 2F(h) + H(h) on the cuts of F(h) 
and H(h), respectively. We have introduced the notation ReF(h) to denote 
the real part on the cut of F(h), and similarly for ReH(h). The principal part 
integral in ( |lj) is thus denoted in ( ^3]) by ReH(h). 

Our object is to now reconstruct the analytic function 2F(h) + H(h) = 
2\nG(h) — H(h) from its behaviour on the cuts. To do this we first need to 
understand the complete structure of cuts. We already know the structure of 
cuts of HQi); it has a logarithmic cut running from h = to h = b, correspond- 
ing to the portion of the density which is saturated at its maximal value of 
one, and a cut from & to a corresponding to the "excited" part of the density, 
where the density is less than one. It thus remains for us to understand the 
cut structure of lnG(/i). 

The function G(h) has two cuts on the physical sheet. The first cut, 
running from b to a, corresponds to the cut of e H ^ h \ the second cut, running 
from cut point c to cut point d, corresponds to the cut of e F ^ . To see whether 
In G(h) has any logarithmic cuts, we first notice from ( |2^) that G(h) is non 
zero everywhere in the complex h plane except possibly at infinity. Thus for 
\nG(h) the only finite logarithmic cut points are at h = b, defined to be the 
end of the flat part of the density (this corresponds to the end of the cut of 
e 11 ^), and possibly the cut point c, defined to be at the end of the cut of 
e F ( h '> . The only remaining question is whether this logarithmic cut starting at 
h = b goes off to infinity or terminates at c. For large h we see from ( |2^) that 
In G(h) = lne + C(i), i.e. there is no logarithmic cut at infinity. We conclude 
that the cut structure of the function In G(h) consists of the cuts corresponding 
to e F (h) and e H ( h > connected together by a logarithmic cut, whose cut points 
are b and c. 

We thus understand the behaviour of 2F(h) + H(h) on all of its cuts. 
Standard methods then allow us to generate the full analytic function 2F(h) + 
H(h). The four cut points, a and b defining the cut of H(h), and c and d 
defining the cut of F(h) are fixed later by boundary conditions. 

Without going into the details of the solution let us present at least one 
final result: the density of weights p(h) of the most probable representation 
R{h}. It can be obtained from 2F(h)+H(h) by using the saddle point equation 
2F(h) + ReH(h) = — \nh and the fact that the resolvent ( [l6]) for the Young 
tableau can be written as H(h) = ReH(h) =p inp(h). To obtain the result 
one has to use all the analytical information listed above. After some tedious 
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calculations we obtain: 
p{h) 

where v and u are denned by 



u 

K 



In 



fl*(gg(" - iv),q) 
8 4 (^(u + iv),q) 



! (a-h)(b-d) 
(a-b)(h-d) 



,k), 



and 



k 2 = 



(a-b)(c-d) 



(24) 



(25) 



(26) 



(o-c)(6-d) 

K and K' are the complete elliptic integrals of the first kind with respective 
moduli k and k' = \/l — k 2 . E is the complete elliptic integral of the second 
kind with modulus k, E(v, k') is the incomplete elliptic integral of the second 
kind with argument v and modulus k' and sn, cn and dn are the Jacobi Elliptic 
functions. The nome q is defined by 

q = e~^. (27) 

To fix the constants a, b, c and d, we expand iirp(h) = 2F(h) + H(h) + lnh 
for large h and compare the resulting power series expansion to that obtained 
from inverting ( [l8|) : 



2F(h) + H(h) = 2 hi e 



,2*4 
v c 2 



1 



1 

7^ 



(28) 



The terms of 0(^j) depend on the as yet unknown positive powers of G in 

Afinal boundary condition is fixed by the normalization of the density. 

Using this information we can calculate the physical quantities of interest, 
for example, the derivative of the free energy with respect to the cosmological 
constant. Denoting the free energy by Z(t, A, /?), Z = e 



-N 2 Z 



we have 



(29) 



where (h) — J dh p(h) h. The result of the last integration, aa well as 
the conditions for a, b, c, d are too bulky to be presented here (see El for the 
details). In the next section we will present already the physical consequences 
of it, namely, the universal expression for the free energy near the flattening 
limit. 
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4 Physical Results and Conclusions 

One needs to do quite tedious calculations to extract the physical results from 
oux solution summerized in the equations ( ^4^ 12^) . The details can be found 
in El and here we present only the essential results: 

1. We found the equation for the critical curve of pure gravity in the 
parameter space (A, /3). It did not show any "flattening" phase transition in 
the whole interval < [3 < oo. Pure gravity appears to be the only possible 
critical behaviour for the large scale fluctuations of the metric (the only infrared 
stable critical point) within our model. This should be true for a big class of 
similar lattice models of pure R 2 gravity in 2 dimensions. For example, if 
we would allow only the deficits of angle being 0, ±7r in the vertices of our 
quadrangulation, this physical conclusion should not change. On the other 
hand, the multicritical behaviours can occur for a more complex parameter 
space (more t* , t couplings to vary) and we could expect some new phase 
structures. Sums over graphs with the sign-changing Boltzmann weights might 
well be necessary for this. 

2. Once we realize that nothing interesting happens in our model for finite 
values of f3 we should concentrate our attention on the double limit (3 — > 0, A — > 
A c = 1 - + 0(/3 2 ) with the double limit parameter x = 1 + (A c - A) 
fixed. Note that this is a very natural, dimensionless scaling parameter; A c — A 
is the continuum cosmological constant with dimension of inverse area and 
(3 controls the number of curvature defects per unit area (and thus also has 
dimension of inverse area). 

We arrive at the free energy in the vicinity of the (locally) flat metric: 

Z(t, A, 0) = ^ [x 6 - + H J x{x > I)*/*] . (30) 

This scaling function exhibits two different behaviours in the two opposite 
limits: 

1. x c — > 1 - the pure 2D gravity regime for any finite positive (3 and 
A-A c (/3): 

Z(t,A,/3)~(A c -A) 5 / 2 (31) 
This corresponds to the standard value of the string susceptibility exponent 

Jstr = -1/2. 

This results mean that even for very big but finite i? 2 -coupling /3 _1 we 
will always find on big enough distances the metric fluctuations obeying the 
pure gravity scaling. In some sense one can say, that the R 2 2D gravity does 
not exist as a special phase of the 2D gravity. The R 2 term can be always 
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dropped from the action in the long wave limit. Our results demonstrate this 
nonperturbatively, starting from a solvable lattice model. 

2. In the opposite limit x — > oo, i.e., (3 — > with A fixed, we obtain: 

*fc^ = £(TzV (32) 

This limit corresponds to a sum over metrics flat everywhere except of a 
finite number of curvature defects. One can easily see from the purely com- 



binatorial calculations that ( 32) corresponds to the sum over all flat metrics 
with only four defects with the deficit of angle tt. Expanding ( |3(]) as a power 
series in (3 (the first term of which is given by ( |3l]) we see that the second 
term in the series expansion corresponds to the the sum over all flat metrics 
with five 7r defects plus one —tt defect, the third to six tt defects plus two — tt 
defects, etc. 

This limit does not appear to describe a new phase of smooth metrics. 
It corresponds instead to the statistics of a dilute gas of curvature defects 
introduced in an otherwise flat metric. It has nothing to do with the 2D gravity 
but it is an interesting statistical mechanical model in itself. More than that, 
one can prove that the dominating graphs are not lattice artifacts: in this 
limit the manifold "forgets" that it is built out of squares. It consists of big 
flat patches with sparsely distributed point-like curvature defects introduced 
ono it in such a way that they close the surface into a manifold of spherical 
topology. If one considers this object as a sphere punctured at the points 
where curvature defects occur, one can show by the direct calculations^, for 
the simplest configurations, that the corresponding moduli parameters do not 
fall on the boundaries of the moduli space but are typically in some general 
position inside the fundamental domain. 

It would be nice to solve this model (and may be more general cases 
including the matter fields) by a continuum approach. In the Liouville theory 
picture the curvature defects correspond to coulomb charges floating in a two- 
dimensional parametrization space. This could be an alternative approach to 
2D quantum gravity in general. 

Let us conclude by noting that our analytical approach to the models 
of dually weighted graphs can be generalized to a big class of other matrix 
models which could be of physical, as well as of the mathematical interest. 
For example, one can formulate a generalized two-matrix model with partition 
function: 



c we thank Paul Zinn-Justin for clarifying this question for us 
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Z= J d N XdT YexpNtr(F(XY) + U(X) + V(Y)) (33) 

where X, Y are the NxN hermitean matrices and F, U, V are some arbitrary 
functions. The character expansion allows us again to reduce the number 
of integration (or summation) variables to ~ N thus making it possible to 
consider it by means of the saddle point methods. Another, more general 
solvable k-matrix model of this type is: 

Z = I ( II dN * X i) ^Ntr{F{\{ X,) + £ U 3 {X 3 )) (34) 

The exp NtrF(l\ K j=1 Xj) factor can be also expended in characters and 
one can then integrate over the angular degrees of freedom of the matrices by 
the use of only the simple SU(N)-orthogonality relations. This approach is at 
the heart of the Itzykson-DiFrancesko formula which we extensively used for 
our present model of lattice 2D R 2 gravity. 

A big question left is how to generalize our approach to include of the 
matter fields on our lattice manifold. The corresponding matrix models con- 
taining the R 2 type couplings can be easily formulated but cannot be solved by 
the character expansion methods presented here. This is not surprising, since 
models of interacting spins on the 2D regular lattice, even the integrable ones, 
are very complicated. Nevertheless, we think that our methods could eventu- 
ally be powerful enough to advance in this direction and provide a missing link 
between two branches of mathematical physics: integrable statistical mechani- 
cal models on regular 2D lattices, on the one hand, and on random dynamical 
lattices, on the other hand. In more physical language, it could provide a link 
between the integrable models of interacting fields with and without the 2D 
quantum gravity fluctuations. 
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